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1 Global theory 

1.1 Statement of the problem 

Let us consider the following operator (quantum Hamiltonian) in R d with 
d = 3 

(1.1) H = H A V = ((hD — A) ■ cr) 2 — V(x) 

where A, V are real- valued functions and V e Jzf i n Jz? 4 , A e J^ 1 . Than 
this operator is self-adjoint. We are interested in Tr - H A V (the sum of all 
negative eigenvalues of this operator). Let 

(1.2) E* = inf EM), 

46^(8(0,1)) 

(1.3) E(A) := (Tr H Ay + n^h' 2 J \dA\ 2 dx^j 

with OA = (djAj) a matrix. Semiclassical asymptotics for such objects were 
studied first in [EFS2] as a first step to consider asymptotics of the ground 
state energy for atoms and molecules in the self- generated magnetic fields 
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[ES, EFS1] which was achieved in [EFS3] where Scott correction term was 
recovered. 

Similarly, this paper is the first step to the recovering sharper asymptotics 
of the ground state energy for atoms and molecules in the self-generated 
magnetic fields, hopefully up to and including Dirac and Schwinger correc- 
tions. To do this we improve theorem 1.1 of [EFS2] (see theorems 2.5, 2.7, 

4.1 and 4.2 below). 

The estimate from above is delivered by A = and Weyl formula with 
an error 0(/) _1 ) as V e c € 2 ' x 1 

(1.4) E* < Wey^ + 0(/?- 1 ); 

where 

(1.5) Weyl(r) = J-h' 3 [{V + tA dx, 

6tt z J 

(1.6) Weyl^ / Td T \Ney\{r) = [ Vl dx. 

J-oo 157r J 

Also for estimates o(h~ 2 ) we need to include into Weyl x the corresponding 
boundary term. The purpose of this paper is to provide an estimate from 
below 

(1.7) E* > Wey^ - O^ 1 ); 

We will use also Weyl(x, r) and Weyl-^x) defined the same way albeit without 
integration with respect to x. 

In this version (v3) we fix several significant errors, include asymptotics 
with o(/? _1 ) remainder and provide a bit more details. 

1.2 Preliminary 

Let us estimate from below. First we need the following really simple 
theorem (cf. 3.1 [EFS1]) 

1 Which means that the second derivatives of V are continuous with the continuity 
modulus 0(| log |x — y|| _1 ), see section 4.5 of [12]). If there is a boundary it does not 
pose any problem as it is in the classically forbidden region. 
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Theorem 1.1. Let V e Jzfi Hi? 4 . Then 

(1.8) E*>-C/T 3 
and either 

(1.9) 4? / \ dA \ 2 dx < Ch ~ 2 - 

Kh 2 J 

or E(A) > c/T 3 . 

Proof. Using the Magnetic Lieb-Thirring inequality (5) of [LLS]) 



;i.lO) / trei(x,x,r) dx > 

-C/T 3 J Vldx-Ch 2 J (V 2 j \dA\ 2 dxy(h-* J V*dx 



\dA\ 2 dx 



we conclude that for any 5 > 

(1.11) E{A) > - C/T 3 - C5 3 h- 3 + (kT 1 - r 1 )/!" 1 J 

which implies both statements of the theorem. □ 
Theorem 1.2. Let V + E Jzfs PI =^ 4 ; « < c/r 1 and 

(1.12) V < -K-\l + \x\) s + K. 
Then there exists a minimizer A. 

Proof. Consider a minimizing sequence Aj . Without any loss of the generality 
one can assume that Aj — y A^ weakly in Jt? 1 and in Jzf 6 and strongly in 
££\ QQ with any p < 6. Then /4 M is a minimizer. 

Really, due to (1.10) non-positive spectra of Ha y are discrete and the 
number of non-positive eigenvalues is bounded by A//,. Without any loss of 
the generality one can assume that Xj^ have limits (we go to the subsequence 
if needed) which are either < or = 0. Here Ay^ are ordered eigenvalues of 
Haj.v- 

We claim that those limits are also eigenvalues and if Xj^, ... , Xj^+r-i 
have the same limit A < 0, it is eigenvalue of at least multiplicity r. Indeed, 
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let Uj t k be corresponding eigenf unctions, orthonormal in J2? 2 . Then in virtue 
of Aj being bounded in J?f 6 and V e Jz? 4 we can estimate 



\\Du hk \\ < K\\uj, k \\\-° ■ \\u hk \\ a < K\\Du hk f- a ■ \\u hk no 

with a > which implies || Duy^H < K. Also assumption (1.12) implies 
that ||(1 + |x|) CT / 2 iv / - i k|| are bounded and therefore without any loss of the 
generality one can assume that Uj tk converge strongly. 
Then 

Mm Tr" H Aj ,v > Tr H Aoo ,v, lim inf J \dAj\ 2 dx > J IdA^] 2 dx 

and therefore E(/4oo) < E* (and then it is a minimizer and there are equalities 
and in particular there no other eigenvalues of Ha^.v- □ 

Remark 1.3. We don't know if the minimizer is unique. Also we do not 
impose here any restrictions on K (which may depend on h) in (1.12) or 
k > 0. From now on until further notice let A be a minimizer. 

Proposition 1.4. Let A be a minimizer. Then 
(1-13) ^AAj{x) = *j := 



- ^2(vjO k (hD k - A k ) x + <s k Oj(hD k - A k ) y )e(x, y,r)\ 



where A = (Ai, A 2 , A3), cr = (oi, o 2 , o 3 ) and e(x,y, r) is the Schwartz kernel 
of the spectral projector of Ha,v- 

Proof. Consider variation 5 A of A and variation of Tr~ (/-/). Note that the 
spectral projector of H is 

(1.14) 9(r - H) = -L r Res K (r - H)" 1 

and therefore 

5Tr6(r-H) = -^ / ReS]jTr(r — H) _1 (6H)(r — H) -1 = 
2vr/ ./ x 

/ T Res K Tr(6H)(r - H)' 2 = -9 T -L f Res K Tr(5H)(r - H)" 1 = 

— d T Tr(5H)0(r — H). 
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Plugging it into 

rd T Tr 0(r - H) = - / Tr d{r - H) dr 

■oo •/ — oo 

and integrating with respect to r we arrive after simple calculations to 

(1.16) 6 Tr~(H) = Tr(5H)0(r — H) = — J 4>(x)8A{x) dx 
where <t>(x) is the right-hand expression of (1.13). Therefore 

(1.17) bE{A) = J (-<D(x) - ^ 2 AA(x))bA(x) dx 

which implies (1.13). □ 
Proposition 1.5. If for re = re* 

(1.18) E* > Weyli - CM 
witt M > C/r 1 /j/ien /or re < re*(l — e ) 

(1.19) 4o [ \8A\ 2 dx< CiM. 

re/? 2 ./ 

Proof. Proof is obvious based also on the upper estimate E* < Wey^ + 
C/r 1 . " □ 

1.3 Estimates. I 

Proposition 1.6. Let (1.19) be fulfilled and let 

(1.20) q = re/M/?i < c 

T/ien as t < c 

(%) Operator norm in ££ 2 of (hD) k d(r — H) does not exceed C for k = 0, 1, 2; 

fnj Operator norm in Jzf 2 of (hD) k ({hD — A) • cr)d(r — H) does not exceed 
C fork = 0,1. 



Proof, (i) Let u = 9(r - H)f. Then ||u|| < ||f || and as 

(1.21) \\A\\ 6 <C\\dA\\<C(KM)h 
we conclude that 

\\hDu\\ < \\{hD - A)u\\ + \\Au\\ < \\{hD - A)u\\ + C||/4|| 6 ||u|| 3 < 
\\(hD - A)u\\ + C{KM)h\\u\\ 3 / 4 \\u\\l //l < 
\\{hD-A)u\\ + C^Hiyll^H/jDivll 1 / 4 < 

||(/,D->\)ii|| + ^||/iDii|| + C?||ii|| 

so due to (1.20) 

(1.22) \\hDu\\ < 2\\{hD - A)u\\ + C\\u\\ 
On the other hand, for B = V x A and r < c 

\\{hD-A)uf < C\\u\\ 2 + {h\B\u, u) < C||u|| 2 + /i||S||||u||2 < 

C\\u\\ 2 + C{KM)h 2 \\u\\ ■ \\u\\ 6 < C\\u\\ 2 + C{KM)h\\u\\\\hDu\\ 

and due to (1.22) we conclude that 

(1.23) \\hDu\\ + \\{hD-A)u\\ < C\\u\\. 

So, for j = 0, 1 statement (i) is proven. 

Further, as h 2 D 2 = (hD - A) 2 + A(hD - A) + AhD - h[D, A] we in the 

same way as before (and using (1.23) conclude that 

\\h 2 D 2 u\\ < C\\u\\ 2 + \\\hD{hD - A)u\\ + \\\h 2 D 2 u\\ 
and therefore 

\\h 2 D 2 u\\ < C\\u\\ 2 + C\\AhDu\\ 

and repeating the same arguments we get ||/j 2 D 2 u|| < C||u||; so for j = 2 
statement (i) is proven. 

(ii) Statement (ii) is proven in the same way. □ 
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Corollary 1.7. Let (1.19) and (1.20) be fulfilled. Then as t < c 

(1.24) e(x,x,r) < Or 3 , | ((/?D - 4) • cr)e(x, y, r)| x=y | < Or 3 . 

Proof. Due to proposition 1.6 operator norms from J?f 2 to ^ of both Q(t — H) 
and ((/?£> — -4) • cr)0(r — H) do not exceed C and the same is true for an 
adjoint operator which imply both claims. □ 

Corollary 1.8. Let (1.19) and (1.20) be fulfilled and A be a minimizer. 
Then for arbitrarily small exponent 5 > 

(1.25) \\dA\\ v i- S < CkH' 1 , 

(1.26) Halloo < /i"H 

where is the scale of Holder spaces and 5 > is arbitrarily small. 

Proof. Really, due to (1.13) for a minimizer HA^H^ < C/t/? -1 . Also we 
know that \\dA\\ < C(kM/) 2 )5 < Ch\ due to (1.20). Then (1.25) holds due 
to the standard properties of the elliptic equations 2 . 

Therefore if at some point y we have |cV\(y)| > /i, it is true in its 
e(/i/?/t _1 ) 1_<5 -vicinity (provided fj, < nh" 1 ) and then 

II^H 2 ^^ 2 ^- 1 ) 3 ^) 

and we conclude that 

^hK' 1 ) 3 ^ < Cnh 2 M ^ n 5 ~ 3S < Cn 4 - 35 h- 1+35 M 

and one can see easily that (1.26) holds due to (1.20) and h~ 1 < M < h~ 3 . 

On the other hand, if fi > nh -1 then we need to take e-vicinity and then 
/i 2 < CnMh 2 < Cti2 where we used (1.20) again. Then (1.26) is proven. □ 

Remark 1.9. (i) It is not clear if it is possible to generalize this theory to 
arbitrary d > 2 with magnetic field energy is given by 

(1.27) ^|(|cVl| 2 -|V-/l| 2 )c/x 

Surely one should use generalized Pauli matrices cr 7 - in the definition of the 
operator. Especially problematic are d > 5. 

2 Actually we can slightly improve this statement. 
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(ii) Therefore while arguments of section 2 below remain valid for d ^ 3, so 
far they remain conditional (if a minimizer exists and satisfies some crude 
estimates) . 



2 Microlocal analysis unleashed 
2.1 Sharp estimates 

Now we can unleash the full power of microlocal analysis but we need to 
extend it to our framework. It follows by induction from (1.25)-(1.26) and 
the arguments we used to derive these estimates that 



so A is "smooth" in e = h scale while for rough microlocal analysis as in 
[BI] and section 2.3 of [12] one needs e = Ch\ log h\ at least. We consider in 
this section arbitrary d > 2; see however remark 1.9. 
Note that 

(2.2) For a commutator of a pseudo-differential operator with a smooth 
symbol and ^ +1 function A a usual commutator formula holds modulo 
0(/? e ' +1 ||(9/4||< (! fe) for any non-integer 9 > 0. 

Proposition 2.1. Assume that 



Then for \a\ < 2, \(3\ < 2, 9 > 1 

(2.5) \F t ^ h - wX T(t)((hD x nhD y fU(x,y,t))\ x=y \ < 

Ch i-d+s T -s + ch- d+e T 2 \\dAy e 

where \ e ^°°([-l, -\) U [\, 1]), Xr(t) = x(t/T), Ch < T < e, |r| < e and 
U(x,y, t) is the Schwartz kernel of e lh ltA . 



(2.1) 



\\dA\\^n-S 



< CnKh- 1 -", 




fi := ||£M||oo.B(x.l) < C , 



V\ > e . 



S 



Proof. Consider first T x 1. First, using standard propagation arguments 
as in [BI] and section 2.3 of [12] one can prove the finite propagation speed 
with respect to x, namely that 

(2.6) \F t ^ h -i TXT {t)(hD x ) a {hD y yU\ < Ch s as \x - y\ > C T, r < c 

with an arbitrarily large exponent s. Further, using (2.6) we prove the finite 
propagation speed with respect to £: 

(2.7) \F^ h -i rX T(t)(hD x ) a (hD y f Vl (hD x ) V2 (hD y )U\ < Ch s + Ch 6 \\dA\\ v e 

as ipi, (f2 G ^o°, dist(supp <p\, supp tp\) > C T, r < c Q 

where the last term in the right-hand expression is due to the non-smoothness 
of A. 

Furthermore, using (2.6) and (2.7) we prove that as (2.4) is fulfilled there 
is a propagation: 

(2.8) \F t ^ h -i TXT {t){hD x ) a {hD y fU\ < Ch s + Ch e \\dA\\ v e 

as |x — y| < €qT, |t| < e, T < e 

which implies (2.5) under additional assumption T x 1; here the last term 
in the right-hand expression is inherited from (2.7). 

Finally, for Ch < T < e we use rescaling t*->t/T,xi->x/T, h i— )■ h/T 
and we use "homogenized" ff® norms. Detailed proof for will be published 
later. □ 

Corollary 2.2. In the framework of proposition 2.1 as \a\ < 2, \/3\ < 2 



(2.9) \F t ^ h - w [Mt)((hD x r(hD y yU(x,y,t))]\ x=y \ < 

Ch l - d + CT 2 h- d+0 \\dA\\v 



where x G ^ , °°([— 1, 1]) and 

(2.10) | [((hD x - A(x)) ■ <r) a {(hD y - A(y)) ■ <rf e(x,y, t)] | - 



WeyU(x)| < Ch'- d + Ch- d+1 ^\\dA\i 
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where 

(2.11) Weyl Q ^(x) = const h~ d [ ((£ - A{x)) ■ <r) a+P d£ 

is the corresponding Weyl expression and 

(2.12) h(x,0={Z-A(x)-ct) 2 -V(x); 

in particular Weyl is as \a\ + \0\ = 1. 

Proof. Obviously summation of (2.5) over C h < \t\ < T and a trivial 
estimate by Ch 1 ~ d of the contribution of the interval |r| < CqIi implies (2.9). 

Then the standard Tauberian arguments and (2.9) imply that (2.10) 
would be correct if we used Tauberian expression with T = T* instead of 
Weyl where we pick up 

(2.13) 7*=emin(l,/)-^ e - 1 )||^|| c ;!). 

Meanwhile the Tauberian formula and (2.5) imply that the contribution of 
an interval {r : \ t\ x 7"} with h < T < T* to the Tauberian expression does 
not exceed the right-hand expression of (2.5) divided by T, i.e. 

Ch i-d+s T -s-i + CTh- d+9 \\dA\ye; 

summation over T* := h 1 ^ 5 < T < T* results in the right-hand expression 
of (2.10). 

So, we need to calculate only the contribution of {t : \ t\ < T*} but one 
can see easily that modulo indicated error it coincides with Weyl a (3 . □ 

Corollary 2.3. (i) If assumption (2.3) is replaced by 

(2.14) /i := Halloo < OT 1+CT 
with a > then 

(2.15) | [((hD x - A(x)) ■ <r) a ((hD y - A(y)) ■ tr) P e{x,y, t)} \ x= - 

Weyl Qi/J (x)| < C~^- d + Cfl-h- d+1 ^\\dA\\l 

and 

(2.16) T* =emm{jl'\]lh-^ e -^\\dA\\:)) 
with ft = max(yU, 1); 
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(ii) As d > 3 one can skip assumption (2.4). 

Proof. As d > 3 (ii) is proven by the standard rescaling technique: x i— > x£, 

3 1 2 

h ^ h = h£~2, A i — ^ /4£~2 with £ = max(e| V\, hi) (see [II] and chapter 5 of 
[12]. 

As n > 1 (i) is proven by the standard rescaling technique x i— >• /ix, 
h h= hfi. □ 

Proposition 2.4. Let k < c, (2.14) be fulfilled, and let A be a minimizer. 
As d = 2 let (2.4) be also fulfilled. Then as 6 e (1,2) 

(2.17) + h^WdAW^e <Ck+ C\\dA\\' 
with 

(2.18) \\dA\\':= sup \\dA\\^ {B{y ^y 

y 

Proof. Consider expression for AA According to (1.13) and (2.15) we get 
then for any 9 G (1,2) 

(2.19) \AA\ + \hdAA\ < Cn{ji + JT^h^ e -^\\dA\\l e ), 
which implies that for any 9' G (1,2) 

h e '- l \\dA\\^ < Cn(ji + ii-h*V- 1 )\\dA\\l 9 ) + C/i, 
and picking up 9' = 9 we conclude that 

(2.20) h^WdAW^e < C(k/2 + K 2 jl- 1 ) + C/i. 

As k < c the right-hand expression does not exceed C(k + /i); then the right- 
hand expression in (2.19) also does not exceed C(n + /i); then ||<9/4||^e-i < 
C(k + /i) and then fi < Ck + C||<9/4||', which implies (2.17). □ 

Having this strong estimate to A allows us to prove 
Theorem 2.5. Let k < c, (2.14) be fulfilled, and let d = 3. Then 

(2.21) E* = Weyli + Of/T 1 ) 
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and a minimizer A satisfies 

(2.22) \\dA\\ < C/cM 
and 

(2.23) + h^WdAW^e < Cith* + Ck. 

Proof. In virtue of (2.9) and (2.17) the Tauberian error with 7x1 when 
calculating Tr H~ v does not exceed 

(2.24) Cjl 2 h 2 - d + C(k + \\dA\\')h 2 - d . 
We claim that 

(2.25) Weyl error when calculating Tr H~ v also does not exceed (2.24). 
Then 

(2.26) E*(A) > Wey^ - Cfl 2 h 2 - d - C(« + ||<94||> 2 - d + k' 1 h 1 ^ \\dA\\ 2 > 

Weyli - Ch 2 - d + ^h 1 - d \\dA\\ 2 

because Ji < C\\dA\\ + 1 due to (2.17). This implies an estimate of E* from 
below and combining with the estimate E* < E*(0) = Wey^ + Ch 2 ~ d from 
above we arrive to (2.21) and (2.22) and then (2.23) due to (2.17). 

To prove (2.25) let us plug A £ instead of A into ei(x,x,0). Then in 
virtue of rough microlocal analysis contribution of {t : 7"* < |t| < e} with 
T* = h 1 ^ 5 would be negligible and contribution of {t : \ t\ < 7"*} would be 
\Ney\ 1 + 0{h 2 ' d ). 

Let us calculate an error which we made plugging A £ instead of A 
into ei(x,x, 0). Obviously it does not exceed Ch~ d \\A — A^^ and since 
\\A - /yioo < Ce e+1 \\dA\\ v e this error does not exceed Ch e+1 - d - AS \\dAye 
which is marginally worse than what we are looking for. However it is good 
enough to recover a weaker version of (2.21) and (2.22) with an extra factor 
h~ Sl in their right-hand expressions. Then (2.17) implies a bit weaker version 
of (2.23) and in particular that its left-hand expression does not exceed C. 

Knowing this let us consider the two term approximation. With the above 
knowledge one can prove easily that the error in two term approximation 
does not exceed Ch 3 ~ d ~ 5 ' with 5' = 1005. 
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Then the second term in the Tauberian expression is 

(2-27) J {(H AV -H A£ , v )el £) (x,y,0))\ y=x dx. 

where subscript means that we plugged A £ instead of A and superscript 
T means that we consider Tauberian expression with T = T* = e. But 
then the contribution of {t : 7~* < \ t\ < T*} is also negligible and modulo 
Ch e+2 ~ d ~* s \\dA\\<ge we get a Weyl expression. However 

(2.28) (H AV - H AeM ) = -2{i - A £ ) ■ (A - A £ ) + \A- A £ \ 2 

and the first term kills Weyl expression as integrand is odd with respect 
to (£ — A £ ) while the second as one can see easily makes it smaller than 
Ch 3 - d - 5 '. Therefore (2.25) has been proven. □ 

Remark 2.6. (i) For d = 2 we cannot skip (2.4) at the stage we did it for 
d > 3. However results of the next section allow us to cure this problem 
using partition-and-rescaling technique. 

(ii) Actually we have an estimate 

(2.29) \dA(x) - dA(y)\ < Ck\x - y|(| log |x - y|| + 1) + C[i. 
Combining with (2.22) we conclude that 

(2.30) Hc^lU < C^ +1 )/( d+2 )| log/)| d /( d+2 V/( d+2 ) 

2.2 Classical dynamics and sharper 
estimates 

Now we want to improve remainder estimate 0(h 2 ~ d ) to o(h 2 ~ d ). Sure, 
we need to impose condition to the classical dynamical system and as 
\dA\ = 0(h cr ) with o > due to (2.30) it should be dynamical system 
associated with Hamiltonian flow generated by H o v : 

(2.31) The set of periodic points of the dynamical system associated with 
Hamiltonian flow generated by Hoy has measure on the energy level 0. 
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Recall that on {(x,£) : H o v (x,^) = r} a natural density d\±r = dxd£ : 
dH\H=T is defined. 

The problem is we do not have quantum propagation theory for H AV as 
A is not a "rough" function. However it is rather regular function, almost ^ 2 , 
and {A-A £ ) is rather small: \A-A £ \ < rj = Ch 2 - 35 and \d{A-A £ )\ < C/? 1 " 35 
and therefore we can apply a method of successive approximations with the 
unperturbed operator H Ae y as long as r/T/h < \f i.e. as T < h 1 ' 46 . Here 
we however have no use for such large T and consider T = 0(h~ s ). 

Consider 

(2.32) F t _ h - wX T{t)U{x,y,t), 

and consider terms of successive approximations. Then if we forget about 
microhyperbolicity arguments the first term will be 0(h~ d T), the second 
0(h^- d r]T 2 ) = Oih 1 "^ 5 ') and the error 0(h- 2 - d r] 2 T 3 ) = 0{h 2 - d - s "). 
Therefore as our goal is 0(h 1 ~ d ) we need to consider the first two terms 
only. The first term is the same expression (2.32) with U replaced by U( £ ). 
Consider the second term, it corresponds to 

(2.33) U{ £) = [ih- 1 f e^ h - lH ^(H AM - H AeM )e it,b ~ lH ^ v dt'] 

Jo 

and then 

(2.34) Tr(e //7_lt » = ih' 1 Tr((H AV - H Ae y) e''"" 1 ^ ^t) 
with 

Jo 

Here [S](x,y) denotes the Schwartz kernel of operator 5. We claim that 

(2.35) \F t ^ h - wX T(t)TrU{ £ ^\ < C V T 2 h- d . 

Here in comparison with the trivial estimate we gained factor h. The proof 
of (2.35) can be done easily by the standard rough microlocal analysis 
arguments and we will provide a detailed proof later. 
Therefore we arrive to the estimate 

(2.36) |F t ^,-i T xr(0Tr((e /t '^ 1 ^^ - e''^ 1 "^) V>) | < Ch 1 ^ . 
On the other hand traditional methods imply that as d > 3 

(2.37) |F t ^- lrXr (t)Tr(e /t '^ 1 ^ ^)| < Ch 1 '" T\i{n T , p ) + C^h 1 -^ 
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where Uj is the set of points on energy level 0, periodic with periods not 
exceeding T, \~\r,p is its p- vicinity, p > is arbitrarily small. 

Theorem 2.7. Let k < c, (2.14) be fulfilled, and let d = 3. Furthermore, 
let condition (2.31) be fulfilled (i.e. [io(\~\oo) = 0)- Then 

(2.38) E* = Weyl* + o^ 1 ) 

where 

(2.39) \Ney\* 1 = \Ney\ 1 + Kh- 1 J vfAV dx 

calculated in the standard way for H o v and a rainiraizer A satisfies similarly 
improved versions of (2.22) and (2.23). 

Remark 2.8. (i) Under stronger assumptions to the Hamiltonian flow one 
can recover better estimates like 0(h 2 ~ d \ \ogh\~ 2 ) or even 0(h 2+s ~ d ) (like 
in subsubsection 4.4.4.3 of [12]). 

(ii) We leave to the reader to calculate the numerical constants here and 
in (3.4), x = K\ — |x2. 



3 Local theory 
3.1 Localization 

The results of the previous section have two shortcomings: first, they impose 
the excessive requirement to k; second, they are not local. However curing 
the second shortcoming we make the way to addressing the first one as well 
using the partition and rescaling technique. 

We localize the first term in E(/\) by using the same localization as in 
[EFS1]: namely we take Tr^Htp) where ip 6 ^(B(0, §)), < ip < 1 and 
some other conditions will be imposed to it later. Note that 

(3.1) Tr(ipHip) > J e 1 {x,x,0)ij 2 {x)dx. 

Really, operator H = HB(-H) + H(l - Q{-H)) where 9(r - H) is a spectral 
projector of H and therefore in the operator sense H > H := HE(Q) and 
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ipHip > ipH if) and therefore all negative eigenvalues of ipHif) are greater 
than or equal to eigenvalues of the negative operator ipH~ip and then 



which is exactly the right-hand expression of (3.1). 

Remark 3.1. (i) The right-hand expression of (3.1) is an another way to 
localize operator trace. Each approach has its own advantages. In particular, 
no need to localize A (see (ii)) and the fact that proposition 1.5 obviously 
remains true (due to corollary 2.3) are advantages of Tr~ (ipHif) )-localization. 

(ii) As Tr~(ipHip) does not depend on A outside of 6(0, |) we may assume 
that A = outside of 8(0, 1). Really, we can always subtract a constant 
from A without affecting traces and also cut-off A outside of 6(0, 1) in a 
way such that A' = A in 6(0, |) and \\dA'\\ < c||cM||g(o,i); the price is to 
multiply k by c _1 - as long as principal parts of asymptotics coincide. 

(iii) Additivity rather than sub-additivity (4.2) and the trivial estimate 
from the above are advantages of Tr^ /-/"^-localization. It may happen that 
the latter definition is more useful in applications to theory of heavy atoms 
and molecules and we will need to recover our results under it. 

Let us estimate from the above: 

Proposition 3.2. Let £(x) be a scaling function 3 and ip be a function such 
that \d a ip\ < c^ _2|a| for all a : \a\ < 2 and < c£ 3i . 
Then, as A = 0, 




(3.2) 




and under assumption (2.31) 



(3.3) Tr-(^H^) = / Weyl*(x)^ 2 (x) dx + o{h~ l ) 




with 



(3.4) Weyl*(x) = Weyl 1 (x) + x 1 /)^ 1 \/]A\/ + x 2 ^ 1 \/]|V\/ 



2 



calculated in the standard way for H o v . 



3 I.e. £> and \d£\ < \. 



4 Such compactly supported functions obviously exist. 
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Proof. Let us consider H = ipHip as a Hamiltonian and let e(x,y, r) be the 
Schwartz kernel of its spectral projector. Then 

(3.5) Tr(^) = J ~ ei {x,x,0)dx = J2 J h{x,x,Q)^dx 

j 

where ip 2 form a partition of unity in M 3 and we need to calculate the right 
hand expression. 

(i) Consider first an e£-admissible partition of unity in 8(0, 1). Let us 
consider 7-scale in such element where 7 = max(d? 2 , h) and we will use 1 
scale in £. Then after rescaling semiclassical parameter rescales 
h 1 — y /7 new = h/'-f and the contribution of each 7-element to a semiclassical re- 
mainder does not exceed Cg(h/'y)~ 1 with g < £ A having the same magnitude 
over element as 7 > 2h. Then contribution of £ element to a semiclassical 
error does not exceed Cg^/j)' 1 x £ 3 *y~ 3 x Ch~ l g£ z £~ A < Chr 1 ^ . 

As £ 2 x h the same arguments work with I replaced by hi and 7 = h and 
effective semiclassical parameter 1. 

Therefore the total contribution of all partition elements in 8(0, 1) to a 
semiclassical error does not exceed Ch -1 . 

(ii) However we need to consider contribution of the rest of IR 3 . Here we use 
7 = ||x|, 1-scale with respect to £ and take g = try~ 2 ; then contribution of 
7-element to a semiclassical error does not exceed Cg^/^y 1 < C^" 1 and 
summation over partition results in C. Thus (3.2) is proven. 

(iii) Note that contribution of zone ||{x : £ < 77} to the remainder does 
not exceed C£ 2 /j _1 ; applying in zone ||{x : £ > i]} sharp asymptotics under 
assumption (2.31) we prove (3.3). □ 

Corollary 3.3. In the framework of proposition 3.2(i), (ii) 

(3.6) E; := inf E^A) < Weylj + Or 1 

and 

(3.7) E; < Weyl* + C/T 1 

respectively with 

(3.8) E^A) := Tr-(V^) + ^ J \dA\ 2 dx. 
Really, we just pick A = 0. 
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3.2 Estimate from below 

Now let us estimate E^(A) from below. We already know that 

(3.9) E^A) > [ e 1 {x,x,0)4> 2 dx+- [ - [ \dA\ 2 dx. 

J Kh 7s(o,i) 

However we need an equation for an optimizer and it would be easier for us 
to deal with even lesser expression involving r-regularization. Let us rewrite 
the first term in the form 

(p(r/L)rd T e(x,x,T)+ / (1 - <^(t/L))t d T e(x, x, r) > 

■oo J — oo 

/ (<p{t/L){t - L) d T e{x, x, r) + (1 - (p{r/L))r d T e{x, x, r)) 

J —CO 

where (p G ^ °°([— 1, 1]) equals 1 in [— j, |] and let us estimate from below 



(3.10) E'^A) =J\J (p{r/L){T - L)d T e{x, x, r)(x)+ 

(1 - <p{t/L)){t - L) d T e{x, x, r))^(x) dx+ 



1 

w 



\dA\ 2 (g(x) + K- 1 )dx. 

6(0,1) 



Proposition 3.4. Let A be a minimizer ofE'^(A). Then 

(3-11) ^2 AA Ax) = <*>J ■= 

((jj(j k (hD k - A k ) x + G k Gj(hD k - A k ) y ) x 

/ \<p{t/L){t - L) Res K (r - H)"V(r - 

J -co 

(1 - £(t//.))t(t - L) Res R (r - H)"V(r - H)" 1 ] (x,y) t/r 

where again [S](x, y) is £/ie Schwartz kernel of 5. 

Proof. Follows immediately from the proof of proposition 1.4. □ 
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y=x 



Proposition 3.5. Let (1.19) and (1.20) be fulfilled. Then as r < c 

(i) Operator norm in J?? 2 of (hD) k (r — H)^ 1 does not exceed C\ lmr| _1 for 
/c = 0,l,2; 

(ii) Operator norm in J2? 2 of (hD) 2 [(hD — A) ■ cr)(r — H)^ 1 does not exceed 
Cllmr]- 1 for k = 0,1, 2. 

Proof. Proof follows the same scheme as the proof of proposition 1.6. □ 
Proposition 3.6. Let (1.19) and (1.20) be fulfilled. Then |<D(x)| < Or 3 . 
Proof. Let us estimate 

(x,y) </t| 

where L < c and y? G ^ °°([— 1, 1]) and a similar expression with a factor 
(t— /_) instead of r; here either T = /, or T = (hD k —A k ) x or T = (hD k — A k ) y . 

Proposition 3.5 implies that the Schwartz kernel of the integrand does 
not exceed Ch~ 3 \ lmr|~ 2 and therefore expression (3.12) does not exceed 
CI 2 x h- 3 L- 2 = Ch- 3 . 

Then what comes out in from the term with the factor (f)(r/h) does 
not exceed Ch~ 3 . 

Representing (1 — 0(r//7)) as a sum of ip(r/L) with L = 2"h with n = 
0, ... , [I log h\\ + c we estimate the output of each term by Ch~ 3 and thus 
the whole sum by Ch~ 3 \ log h\. 

To get rid off the logarithmic factor we rewrite (r — H)~ l ip(r — H)^ 1 as 
— d(r — H)~ x il) + (r — H)~ 2 [h, ip](r — H) -1 ; if we plug only the second part 
we recover a factor h/L where h comes from the commutator and 1/L from 
the increased singularity; an extra operator factor in the commutator is not 
essential. Then summation over partition results in Ch~ 3 . 

Plugging only the first part we do not use the above decomposition but 
an equality Res R (r — H)^ 1 dr = d T d(r — H). □ 

Corollary 3.7. Let (1.19) and (1.20) be fulfilled and A be a minimizer. 
Then (1.25) and (1.26) hold. 

Proof. Proof follows the proof of corollary 1.8. □ 
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(3.12) 



Mt/L) Res K T(r - H)"V(r - H)~ 



Now we can recover both proposition 2.4 and finally theorems 2.5 and 2.7: 
Theorem 3.8. Let (1.19) and k<c be fulfilled. Then 

(i) The following estimate holds: 

(3.13) E; = J WeyUM^x) dx = O^ 1 ) 

and and a minimizer A satisfies (2.22) and (2.23); 

(ii) Furthermore, let assumption (2.31) be fulfilled (i.e. (^(rioc,) = 0). Then 

(3.14) E; - J WeylKx)^ 2 (x) cfr = o^ 1 ) 

and a minimizer A satisfies similarly improved versions of (2.22) and (2.23). 



4 Rescaling 
4.1 Case ac < 1 

We already have an upper estimate: corollary 3.3. Let us prove a lower 
estimate. Consider an error 

(4.1) (| Weyl^x^c/x-E^)). 
Obviously Tr~ is sub-additive 

(4.2) Tr-(5>^y) > E Tr "(^ H ^) 

and therefore so is Ey,(/4). Then we need to consider each partition clement 
and use a lower estimate for it. While considering partition we use so called 
ISM identity 

(4.3) E^ = 1 =► W = E(^ + ^ H '^^])- 
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In virtue of theorem 1.1, from the very beginning we need to consider 

(4.4) M = K p h~i' a 

with a = §, P = and k < c. But we need to satisfy precondition (1.20) 
which is then 

(4.5) K^ +1 h~ a < c. 

If (4.5) is fulfilled with a = we conclude that the final error is indeed 
0(/? _1 ) or even o(/? _1 ) without any precondition. 

Let precondition (4.5) fail. Let us use 7-admissible partition of unity ipj 
with ipj satisfying after rescaling assumptions of proposition 3.2. 

Note that rescaling x y x/7 results in h 1— > /j new = /j/7 and after 
rescaling in the new coordinates ||<9/4|| 2 acquires factor 7 rf ~ 2 and thus factor 
K~ 1 h 1 ~ d becomes Kr 1 h 1 ~ d ^ d ~ 2 = K~^h^ with /c (->• /t new = ^7. 

Then after rescaling precondition (4.5) is satisfied provided before rescal- 
ing ^+i/7-ay*+/3+i < c. Let us pick up 7 = /,«/("+£+!). 

Obviously if before rescaling condition (4.5) failed, then /) C 7 < 1. 

But then expression (4.1) with ^ replaced by ipj does not exceed C(h/ / y) 2 ~ d 
and the total expression (4.1) does not exceed C/? _1 7~ 2 = Ck 13 ' h~^~ a ' with 

= 2(/3 + l)/(a + /3 + 1), «' = _! + 2a/(o! + /5 + 1); 

So, actually we can pick up M with a, (5 replaced by a', (3' and we 
have a precondition (4.5) with these new a', (3' and we do not need an old 
precondition. Repeating the rescaling procedure again we derive a proper 
estimate with again weaker precondition etc. 

One can see easily that a' + j3' + 1 = | and therefore on each step 
a + p + 1 = I and we have recurrent relation for a 1 : a' = — \ + |«; and 
therefore we have sequence for a which decays and then becomes negative. 
Precondition (4.5) has been removed completely and estimate M = 0(h^ 1 ) 
has been established. After this under assumption (2.31) we can prove even 
sharper asymptotics. Thus we arrive to 

Theorem 4.1. Let d = 3, V £ < ^ 2,1 , k < c and let ip satisfy assumption of 
proposition 3.2. Then 

(i) Asymptotics (3.13) holds; 
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(ii) Further, if assumption (2.31) is fulfilled then asymptotics (3.13) holds; 

(Hi) If (3.13) or (3.14) holds for E^(A) (we need only an estimate from 
below) then \\dA\\ = 0(kH)2) or \\dA\\ = o(kH)2) respectively. 

4.2 Case 1 < k < h~ l 

We can consider even the case 1 < k < h" 1 . The simple rescaling-and- 
partition arguments with 7 = lead to the following 

(4.6) As 1 < k < h~ l remainder estimate 0(k 2 /7 _1 ) holds and for a mini- 
mizer \\dA\\ 2 < Cn 3 h. 

However we would like to improve it and, in particular prove that as k 
is moderately large remainder estimate is 0(/7 _1 ) and even o(h~ l ) under 
non-periodicity assumption. 

Theorem 4.2. Let d = 3, V G < ^ 2,1 , and let ip satisfy assumptions of 
proposition 3.2. Then 

(1) As 

(4.7) k < K* h := eh~*\ log 
asymptotics (3.13) holds; 

(ii) Furthermore as k = o(n* h ) and assumption (2.31) is fulfilled then asymp- 
totics (3.14) holds; 

(Hi) As 1 < k < ch^ 1 the following estimate holds: 

(4.8) ' E ^~/ We y'iW^ 2 dx \ ^ Ch~ 3 {Kh)^\ \ogKh\ 2 

Proof, (i) From (2.20) we conclude as k > c that h^ldA^e < C/t(/t + /2). 
Then using arguments of subsection 2.2 one can prove easily that for k < h a ~2 

\F t ^ T XT(t)(hD x r(hD x f(U(x,y,t)-U £ (x,y,t)-U' £ (x,y,t))\ < Ch l - d 

where we use the same 2-term approximation, T = ejr 1 . Let us take then 
x = y, multiply by e~ d ip(e 1 (y — z)) and integrate over y. Using rough 
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microlocal analysis one can prove easily that from both U £ (x, y, t) and U' £ we 
get 0(h 1 ^ d ) and in the end of the day we arrive to the estimate \AA £ \ < Cnjx 
which implies 

(4.9) \d 2 A £ \ < CkJi\ \ogh\ + C/i 

where obviously one can skip the last term. Here we used property of 
the Laplace equation. For our purpose it is much better than |<9M e | < 
Ck 2 \ \ogh\ + Cfi which one could derive easily. 

Again using arguments of subsection 2.2 one can prove easily that 

(4.10) | Tr(VH^) " Tr(^ £ ^)| < Cji 2 h 2 - d 
and therefore 

(4.11) | Tr(^H- u V) - J Weyl 1 (x)^ 2 (x) dx\ < Cjl 2 h 2 - d 

and finally for an optimizer 

(4.12) \\dA\\ 2 < Cntfh. 

Here \i and p, were calculated for A, but it does not really matter as due 
to \d 2 A\ < Cn 2 h- 5 we conclude that \dA - dA £ \ < Ck 2 Ii- 5 e < C due to 
restriction to k. 
Then, as d = 3 

(4.13) fi 2 (fi/nfl\ \ogh\) 3 < Kfl 2 h 

and if /i > 1 we have jl = /i and (4.13) becomes k~ 3 \ log h\~ 3 < Cnh which 
impossible under (4.7). 

So, fi < 1 and (4.13) implies (3.13) and (4.12), (4.13) imply that for an 
optimizer \\dA\\ < C(kH)? and fi < Cn A h\ \ogh\ d . So (i) is proven. 

(ii) Proof of (ii) follows then in virtue of arguments of subsection 2.2. 

(iii) If n* h < k < h^ 1 we apply partition-and-rescaling. So, h i— >■ h' = h/^y and 

4 1 

k i — y k! = K'y and to get into (4.7) we need 7 = e«: — 3 /j — 3 1 \og[nh)\ leading 
to the remainder estimate Ch^ 1 ^ 2 which proves (ii). □ 

Remark 4-3. In versions 1 and 2 (vl and v2) we lost a logarithmic factor in 
(4.7) and (4.8) (it was h~* and h~iKi). 
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